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We consider (2+l)-dimensional topological quantum states which possess edge states described 
by a chiral (l+l)-dimensional Conformal Field Theory (CFT), such as e.g. a general quantum Hall 
state. We demonstrate that for such states the reduced density matrix of a finite spatial region of the 
gapped topological state is a thermal density matrix of the chiral edge state CFT which would appear 
at the spatial boundary of that region. We obtain this result by applying a physical instantaneous 
cut to the gapped system, and by viewing the cutting process as a sudden "quantum quench" into 
a CFT, using the tools of boundary conformal field theory. We thus provide a demonstration of the 
observation made by Li and Haldane about the relationship between the entanglement spectrum 
and the spectrum of a physical edge state. 
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Topological phases of matter are gapped quantum 
states which cannot be adiabatically deformed into a 
completely 'trivial' gapped system such as a trivial 
band insulator, without crossing a quantum phase transi- 
tion. They are not characterized by symmetry breaking, 
but instead by certain global topological properties such 
as the presence of (topologically) protected edge states 
and/or a ground state degeneracy which depends on the 
topology of the surface on which the state resides [l]. 
Topological states of matter of this kind which have 
been discovered in nature include the integer and frac- 
tional quantum Hall states Q, and the recently discovered 
time-reversal invariant topological insulators [" 



Quantum entanglement is a purely quantum mechan- 
ical phenomenon which has no classical analog. For any 
pure quantum state (typically the ground state) of a sys- 
tem consisting of two disjoint subsystems A and B, sub- 
system A can be described by a density matrix obtained 
by 'tracing out' the degrees of freedom in B. This den- 
sity matrix provides complete information about the en- 
tanglement properties of the initial pure state between 
the two subsystems. Quantum entanglement provides 
an alternative characterization of the properties of the 
many-body system, in particular for topological states 
of matter which cannot be described by conventional 
probes such as order parameters, correlation functions, 
etc.. Q For example, as discovered by Levin and Wen, 
and by Kitaev and Preskill0, Q, the entanglement en- 
tropy of a topologically ordered state in a region of lin- 
ear size I in two-dimensional position space contains a 
universal ^-independent term, the 'topological entangle- 
ment entropy' (TEE), which is a characteristic of the 
topological order of the state. However, the TEE is not 
a complete description of a topological state of matter, 
since distinct topologically ordered states can have the 
same TEE. More complete information about a topolog- 



ical state of matter can be obtained from the eigenvalue 
spectrum of the reduced density matrix, often referred to 
as the entanglement spectrum [9(. In general, the density 
matrix pa describing the entanglement between a subsys- 
tem A and the rest of the system can be written in the 



form of pa 



with He a Hermitian operator. This 



is so far nothing but the definition of He- In general, the 
so-defined 11 entanglement Hamiltonian!^ He is different 
from the physical Hamiltonian of the system. (In this 
form, the operator He appears formally in a way anal- 
ogous to that of the physical Hamiltonian f3H in a sys- 
tem in thermal equilibrium, which has a thermal density 
matrix p oc e~@ H at temperature T = l/fce/3.) One im- 
portant physical feature of the so-defined entanglement 
Hamiltonian He is that its low-energy eigenstates corre- 
spond to those states in subsystem A, appearing in the 
Schmidt-decomposition of the initial pure state, which 
are most entangled with the rest of the system. 

The focus of the present article is a remarkable obser- 
vation made recently by Li and Haldane Q and in subse- 
quent work for topological phases whose physical Hamil- 
tonian possesses low energy states at an open boundary 
('edge states'). This includes fractional quantum Hall 
states non- interacting topological insulators [l3l 
and the Kitaev honeycomb model[l5|]. It was found 
that for those systems the low-energy 'edge' states of the 
physical Hamiltonian at an actual open boundary of sys- 
tem A are in one-to-one correspondence with the low- 
lying eigenstates of entanglement Hamiltonian He {i.e. 
with the most entangled states). However, except for 
systems which can be reduced to non-interacting fermion 
problems [l3, 15, 16|, such a correspondence between en- 
tanglement spectrum and the edge state 1 spectrum of the 
physical Hamiltonian has only been supported by nu- 
merical evidence. No general argument for the vali dity 
of such a correspondence has been presented so far[41(. 
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It is the purpose of the present Letter to demonstrate the 
general validity of this correspondence. 
General setup- In this Letter, we show that for a generic 
(2 + l)-dimensional topological state which possesses 
edges states described by a conformal field theory, the 
entanglement Hamiltonian He is proportional to the 
Hamiltonian Hl of a physical chiral (say L-moving) edge 
state appearing an actual spatial boundary of subsystem 
A, in the long-wavelegnth limit and in any fixed topologi- 
cal sector. For example, our conclusion applies to all the 
Abclian and non-Abelian quantum Hall states described 
by Chern-Simons effective field theories in the bulk[I3- 
[20|. In order to relate the entanglement spectrum and 
spectrum of the physical edge state Hamiltonian, we 
consider a bipartition of the toplogical state on a cylin- 
der into two parts A and B as shown in Fig. Q] (a) . The 
(physical) Hamiltonian H can be written in the form 

H = H a +Hb+H A b (1) 

where Ha and Hb denote the Hamiltonians in discon- 
nected regions A and B, respectively, each of which has 
(two) open boundaries. The term Hab couples regions 
A and B across their joint boundary. For example, 
for a 2D gapped tight-binding model H = Yluj) c l^ij c j 
realizing[2lj the integer quantum Hall effect, the term 
Hab contains all the electron hopping terms across the 
boundary between A and B. 

Now we consider a deformed Hamiltonian containing a 
parameter A E [0, 1] (similar to Ref. (22l|): 

H(X) =H A + H B + XHab (2) 

By construction, H(X = 0) is the Hamiltonian of the 
two decoupled cylinders A and B, and H(X = 1) is the 
Hamiltonian of the whole cylinder A U B. Since we are 
interested in such topological states which possess chiral 
edge states, the Hamiltonian H(X = 0) will have chiral 
and anti-chiral edge states propagating at the boundary 
between regions A and B, as shown in Fig. [T] (b). When 
A ^ 0, the term XHab introduces a coupling between the 
regions A and B. Denote the bulk gap of the Hamilto- 
nian H = H(X = 1) by Eg. When the coupling A is 
small enough such that the energy scale of the coupling 
term XHab is much smaller than the bulk gap E s , the 
gapped bulk states described by Ha and Hb are almost 
entirely unaffected by the coupling term XHab, whose 
main effect is then to induce an inter-edge coupling be- 
tween the chiral and anti-chiral edge states. Since each 
individual edge state is described by a chiral conformal 
field theory (CFT), the theory of the two edges between 
regions A and B is described by a non-chiral conformal 
field theory. Thus at low-energy the coupling term XHab 
between regions A and B is reduced to a local interaction 
in the CFT describing the dynamics of the two coupled 
edges. If this interaction is a relevant perturbation of 
the CFT describing the decoupled edges, then the two 



counterpropagating edges will be gapped for arbitrarily 
small coupling A. Thus we expect that in this case the 
system described by the Hamiltonian H(X = 1) to be 
adiabatically connected to that described by H(X) for a 
small but non- vanishing value of A. In this case, the en- 
tanglement properties of H(X = 1) are expected to be 
qualitatively the same as those of H(X) with a small A. 
The latter describes the entanglement between the left- 
and the right-movers of the edge state CFT. Let us as- 
sume for simplicity that this coupling is relevant in the 
renormalization group (RG) sense. [42j 

Below we will solve this entanglement problem for the 
edge state CFT, by mapping it to a problem of a quan- 
tum quench. We then solve the latter (quantum quench) 
problem in the standard manner by using the work of 
Calabrese and Cardy[23|,[24| which employs the methods 
of boundary conformal field theory (BCFT) [25| ■ 
Reduced density matrix of the edge CFT- Next we study 
the entanglement properties of the Hamiltonian H (A) for 
small values of A, which, as explained above, amounts to 
the study of the (1 + 1) dimensional problem of coupled 
edge states, 

Hodge (A) = Hl + Hn + XH int (3) 

Here, Hl and Hn denote the Hamiltonians of left-moving 
(L) and right-moving (R) edge states, and AH; nt a rel- 
evant inter-edge coupling. The left- and right- moving 
edge states are the low-energy excitations of the subsys- 
tem in regions A and B, respectively. Again, the en- 
tanglement properties between the subsystems A and 
B are reduced to those between left and right mov- 
ing (1 + l)-dimensional edge states. If we denote the 
ground state of the Hamiltonian H e d ge (A) from Eq.Q 
by G, then our goal is to obtain the density matrix of 
the left-moving edge state subsystem defined by 

Pl =Ttb(\G)(G\), (4) 

where Tr# denotes the trace over the right-moving edge 
state degrees of freedom. In general, the ground state 
\G) will depend on all the details of the coupling between 
the right- and the left- moving edges states. However, 
due to the gapless nature of H e dge(A = 0) describing 
the decoupled edges, certain universal properties can be 
inferred without reference to any detailed features of this 
coupling in the long- wavelength limit. 

In order to understand the entanglement properties 
of the state G), we relate them to another problem 
- the "quantum quench" problem. Consider a "quan- 
tum quench" of the system composed of the coupled 
edges, Eq. [3] For all times t < the system is in the 
ground state \G) of the Hamiltonian H e d ge {Xo) with non- 
vanishing coupling Ao 7^ between the edges. At time 
t = the coupling Ao between the edges is suddenly 
switched off, so that A = for t > 0. After the quantum 
quench, the left and right moving edge states evolve inde- 
pendently with the Hamiltonian H cc \ gc (X = 0) = Hl+Hr 
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FIG. 1: (a) A topological state on a cylinder with a bipartition 
into two regions A and B. (b) The deformed system (see 
text) with the coupling between A and B regions weighted 
by a factor A £ [0, 1]. The system can be understood as 
two cylinders A and B, with edge states propagating along 
the boundary between A and B, coupled by an inter-edge 
coupling, (c) For small enough A, the coupling between the 
gapped bulk states can be neglected, and the problem can 
be reduced to an inter-edge coupling problem described by 
a (1 + 1) -dimensional conformal field theory with a relevant 
coupling Aifmt- 



of the decoupled edges. Space- and time-dependent cor- 
relation functions after a sudden quench, as above, have 
been studied extensively by Calabrese and Cardvj23ll24|- 
who applied BCFT to obtain general properties of such 
correlation function in the long-time and long- wavelength 
regime. This is relevant for our purpose because the 
density matrix pl is uniquely determined by the set of 
all the equal-time correlation functions of operators with 
support solely on the left-moving edge, 



C(t,{ Xi }) = 

( G \ e ^ L +H R ) dL i(xi y _ dL n{xn)e -U(H L +H R) |G) s 

i- ltHL p L e ltHL 6 L , 1 (x 1 )...6 L , n (x n )\ , (5) 



(All the coordinates x\, x%, ■■■X n reside entirely on the 
left-moving edge.) In the quantum quench problem, 
the ground state G) of the coupled edge Hamiltonian 
H e dge(^o 0) represents an initial condition at time 
t = for the evolution with the gapless (critical) decou- 
pled edge system Hamiltonian Hl + Hr at subsequent 



2G] 



times t > 0. This initial state can be viewed |23|, |2 
as a boundary condition on the gapless theory of the 
right and left moving edges. It can thus be described 



using the methods of boundary critical phenomena [27 1. 
This can be achieved in any dimension of space by an- 
alytical continuation of the real-time Kcldysh contour to 
imaginary time, and subsequent exchange of the roles of 
space and imaginary time (possible due to the underly- 
ing effective relativistic invariance of the low-energy edge 
state theory). Moreover, in the present case of a one- 
dimensional edge, the resulting boundary condition can 
be analyzed by using the powerful tools of BCFT 25 1. 
The key result that we will use from the theory of bound- 
ary critical phenomena is that an arbitrary boundary 
condition on a gapless bulk theory will always renor- 
malize at long distances into a scale invariant bound- 



ary condition [23j, |28|, |29|. Moreover, in the case of a 
(1+1) conformal bulk theory, such as the one describing 
the one-dimensional edges, any scale invariant boundary 
conditions must be one of a known list of conformally 
invariant boundary conditions j25j|. Consequently, as 
emphasized in Ref. 



23j, 



in the long wavelength limit 
the correlation functions at a general boundary condi- 
tion described by a general state \G) are equal to those 
at a conformally invariant boundary condition described 
by a state |G*) which represents a (boundary) fixed point 
to which the boundary state G) flows under the renor- 
malization group (RG). The difference between \G) and 
|G*) can be represented by an imaginary time evolution 
operator, 



\G) ~ Z - l l 2 e-^ HL+H ^ \G,) 



(6) 



where r > is the so-called extrapolation length[23j, |27j 
and stands for the RG "distance" of the general bound- 
ary state \G) to the conformal boundary state G»). 
Z = (G*| e -2T (H L +H R ) is a noma ii za tion factor. 
Physically, the energy scale 1/to is determined by the 
energy gap E(Xq) induced by the coupling term Ao-ffmt 
between the edges. 

In a so-called rational CFT such as the one under con- 
sideration, all conformal invariant boundary states |G*) 
are known [25| to be finite linear combinations of so-called 
Ishibashi states [3(3] which have the form 

|G*, a ) = 

00 d a (n) 

= /2z2 \k(a,n),j;a) L ®\~k(a,n),j;a) R . (7) 

n=Q j=l 

Here a denotes a topological sector in the underlying 
topological theory, i.e. a topological flux threading the 
cylinder in Fig. ((T|a)), which is represented in the CFT 
describing the edges by a primary state of a correspond- 
ing conformal symmetry algebra (Virasoro or other) of 
conformal weight h a . (a denotes the conjugate sector 
and state of conformal weight h a = h a .) The label 
a runs over all possible particle types of the topological 
state [3lj]. Here k(a,n) = 2ir(h a +n)/l denotes the mo- 
mentum, where where I is the circumference of the edge 
of the cylinder; j = 1, 2, ..d a (n) labels the elements of 
an orthonormal basis in the subspace of fixed momen- 
tum k(a,n). Notice that the left- (right-) moving edge 
system only contains excitations with positive (negative) 
momentum. We note that the state in Eq. [7] is and exam- 
ple of a so-called maximally entangled state. The explicit 
form, Eq. [71 of the Ishibashi states |G* ia ), resulting from 
conformal invariance, is of great help in determining the 
form of the reduced density matrix for the left-moving 
edge. Upon directly combining Eqs. ([6]) with (|7J) one 
obtains 



\Ga) 
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E 



-2r vk(a,n) d ^ n ) 



ZV2 

which yields the following form of the density matrix of 
the left-moving edge upon tracing out the right moving 
edge, 

PLa = Tr R (\G a ){G a \)~ 

00 e -iT vk(a,n) 

E 7 E \ k ( a i n )d\a) L (k{a,n) 1 j\a\ L = 

3=1 

(8) 



Z 

Z- l P n e- iToHL P n 



Here we have used the linear dispersion Hl \k,j;a) L = 
vk(a,n)\k,j;a) L , H R \-k,j;a) L = vk\-k,j;a) R where 
v is the edge state velocity and k stands for k(a,n). 
The label a indicates that pLa is an operator defined in 
the topological sector corresponding to topological flux 
a threading the cylinder, and P a is the projection oper- 
ator onto that sector of the Hilbert space of the CFT. 
In cylinder geometry there is no entanglement between 
different topological sectors (denoted by different labels 
a). 

Eq. ([8]) is the central result of this work, which demon- 
strates that the entanglement between left-moving and 
right-moving edge states in a CFT induced by a rel- 
evant coupling is always characterized by a "thermal" 
density matrix within a fixed topological sector (or pri- 
mary state, in the CFT context). In other words, in 
each topological sector the "entanglement Hamiltonian" 
H E = — log pl = 4:T H L + \ogZ is proportional to the 
Hamiltonian Hl of a physical edge up to a possible shift 
of the ground state energy in that sector which ensures 
the proper normalization of the density matrix as a prob- 
ability distribution. Our result demonstrates not only 
that the excitation energies of the entanglement spec- 
trum are the same as those of the spectrum of the Hamil- 
tonian of the edge state of the topological system appear- 
ing (by assumption) at a physical boundary of region A, 
in the long-wavelength limit modulo a global rescaling, 
but also that the most entangled states are in one-to-one 
correspondence with the low-energy edge states which 
occur at this boundary. 

Example: Free fermions- A simple example in which 
the general notions, developed in the preceeding part 
of this article, can also be illustrated using elemen- 
tary many-body techniques is that of the 2D integer 
quantum Hall (IQH) state (and more generally non- 
interacting topological insulators). This state can be 
described by a free fermion theory, the entanglement 
properties of which have been studied extensively in the 
literature [lU Qjl [32|. However, it is still helpful to present 
the results here as an illustration, in the language of the 
much more general formulation obtained above. The 
edge states of an IQH state with integer filling frac- 
tion v — N consist of N flavors of non-interacting chiral 



fermions. For simplicity, we consider an IQH state with 
filling fraction N = 1, whose edge state dynamics is 
governed by the Hamiltonian 



H L = "£v 

k 



kc\c k , H R = - vkd\d k 



(9) 



The simplest inter-edge coupling term is a single-particle 
inter-edge tunneling 



#hrt = Eg \ C l dk + d l Ck 



(10) 



with E s the bulk gap which acts as a high-energy cut- 
off scale for the edge theory. The coupled Hamiltonian 
Hl + Hn + XH[ nt is a free Fermion Hamiltonian which 
can be diagonalized by a unitary transformation to Hl + 
Hr + H int = Z)jfe, s =±i E kll s lks with the gapful energy 

dispersion E^ = ^Jv 2 k 2 + E 2 . Here 7^ (i = 1,2) are 
quasiparticle annihilation operators. The ground state 
\G) of this gapped system is determined by the conditions 
lk,i\G) — (i — 1,2). One obtains[43| the following 
explicit expression for \G) (unnormalized): 



\G) 



-H. 



\G*) 



(11) 



with 



|G*} = cxpi -J2 {c\d k +d 

I k>0 




and H e ~ 2^~{Hl + Hr) in the long wavelength limit. 

The operators c\dk and d)_ k c-k with k > create quasi- 
particle excitations of the system of the two edges, so 
that |G*} is an equal- weight superposition of all quasi- 
particle excitation states in the massless theory; this is 
nothing but the Ishibashi state for the Free fermion CFT 
(in the sector without topological flux) . Thus, with this 
form of H e , we recover correctly (in the long wavelength 
limit) the general relation ([5]); the extrapolation length 
is t = 1/2E & . As expected, the energy scale l/r is 
determined by the energy gap 2E g of particle-hole exci- 
tations. 

Conclusion and discussion- In conclusion, we have 
demonstrated that for a generic (2 + 1) dimensional topo- 
logical state possessing gapless edge states which are de- 
scribed by a chiral CFT, the reduced density matrix of of 
a region A obtained by tracing out the rest of the system 
always has, in the long wavelength limit, the same form 
as a thermal density matrix of region A with an open 
physical boundary. Besides topological states, our anal- 
ysis also applies to other systems described by coupled 
CFTs. In particular, our result provides an explanation 
of the recent numerical and analytical results on the en- 
tanglement spectrum of coupled spin chains [U and cou- 
pled Luttinger liquids [34|. Since the relationship between 
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a general boundary state and a scale invariant boundary 
condition which is the endpoint of the RG flow also holds 
for higher dimensional scale invariant bulk theories 27J , 
we expect that our result will generalize to higher dimen- 
sional topological states, such as (3+1) dimensional topo- 
logical insulators, and especially the fractional topolog- 
ical insulators 35l438[ which cannot be analysed using 
free fermion methodsjl3l. llH ] . Details of this generaliza- 
tion will be left for future work. 

We would like to note that the reduced density matrix 
flSJ) in the topological sector a yields an entanglement 
entropy of the form S = — Tr (p^ log pi,) = aL — Stopo, 
with Stono = log(-D/d Q ) the topological entanglement 
entropy 0, HJ. Here d a is the quantum dimension of the 
quasi-particle of type a, and D = y$Ja i s the total 
quantum dimension. This relation to the topological en- 
tropy has been noticed in Ref. Q, though in that work 
the form (|SJ) of the density matrix was taken as an as- 
sumption. The present paper proves this assumption. 
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Supplentary Material I: More detailed discussion on irrelevant inter-edge coupling 

In the main text we have mainly studied the cases in which the interaction term A/fint in Eq. ([3]) is a relevant 
coupling in the edge CFT, so that a gap is induced once a finite coupling A is turned on. In this appendix, 
we will discuss the cases in which the interaction term \H{ nt is irrelevant, and provide arguments that our result on 
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entanglement spectrum still holds in this case. 

For concreteness of the discussion, we study the simplest fractional quantum Hall state-1 / m Laughlin state as an 
example. In the cylinder geometry shown in Fig. [T] the two edges in the middle are described by a Luttingcr liquid 
(see e.g. [391 ]) 



TTXi Tft 

C Q = — (d t - vd x ) 4>Ld x cf)L + 7T- (-d t - vd x ) <t>ndx4>R 

ZTT Z7T 



(13) 



The inter-edge interaction can be written as[4C 



A cos 



(14) 



Physically, electron inter-edge tunneling corresponds to R = 1. When a forward scattering term such as gd^^d^^^ 
is turned on together with the back-scattering term, the compactification radius R can deviate from 1. For fractional 
quantum Hall edge states m > 1 is an odd integer, and the coupling A is irrelevant if R = 1, as shown in Fig. [2] If A 
is continuously tuned starting from A = 0, the system remains gapless until A = A c where a phase transition occurs 
and the system becomes gapped. By construction, the physical system (Fig. [T] (a)) without edge between A and B 
region corresponds to some value A > A c in the gapped phase, as shown in Fig. [2] by point A. For A < A c the coupling 
is irrelevant, but the compactification radius R of the theory is renormalized. In other words, for A < A c the long 
wavelength behavior of the coupled system is still described by a CFT but it is different from the original CFT in 
Eq. (|13p . Thus naively it seems that our derivation in the main text does not directly apply to this theory However, if 
we also turn on the forward scattering which increases R, the scaling dimension of A can be tuned to the region where 
A is relevant, as shown in Fig. [5] by point B. At this point Eq. (J5J applies, and the reduced density matrix can be 
written as pi, = Z~ 1 e~ 4:T ° HL in each given topological sector. Now consider a path in the parameter space connecting 
A and B parameterized by (X(t),R(t)) with t E [0,1]. As long as the path stays in the gapped phase, for each 
point on the path the system has a gapped ground state \G(X(t), R(t))), which determines a reduced density matrix 
PL(X(t),R(t)) = Ttr \G(X(t), R(t))) (G(X(t), R(t))\ correspondingly. Since no phase transition occurs along the path, 
the wavefunction of the state \G(X(t), R(t))) is a smooth function of A and R, so that pi,(X(t), R(t)) is also smooth. 
Consequently, there is a one-to-one correspondence between the states in the low- lying entanglement spectrum of pl 
at the two points A and B. In other words, the low-lying entanglement spectrum of point A contains the same states 
as the chiral CFT of the edge theory. The eigenvalues of He = — log pL(X(t), R(t)) can change continuously during 
the deformation, but the long-wavelength limit of the dispersion relation has to remain linear since the deformation 
is smooth. Thus in the long wavelength limit, the reduced density matrix p^ at A point still has the form of Eq. ([8]), 
with generically a different tq from B point. 




FIG. 2: Illustration of the RG flow of model (|13[) and (|14[) with two parameters A and R. Any two points A and B in the 
gapped phase can be connected by a continuous path (red dash line). By this continuous deformation one can show that the 
entanglement spectrum of system A is qualitatively the same as system B. The entanglement Hamiltonian of system B can be 
obtained using the approach in the main text since A is relevant. 



The argument above applies to more generic CFTs as long as such path in the parameter space of the theory exists, 
which smoothly connects the parameter region where the coupling is irrelevant to the region where the coupling is 
relevant. In other words, when the coupling XH m t is irrelevant, our conclusion on the relation between entanglement 
spectrum and edge state spectrum still holds as long as there is a marginal coupling in the CFT which can tune the 
scaling dimension of H mi continuously to make it relevant. 
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Supplementary Material II: Detailed derivation of the free Fermion density matrix 

The free fermion Hamiltonian given by Eq. (|9|) and (fTU|) can be diagonalized in the following form: 

H L + H R + XH int = EklLlks (15) 



fe,s=±l 



with E k = Jv 2 k 2 + E%, 



I E k +vk 

lk+ \ I U k -Vk \ I Cfe \ I U k \ (A/ 2E k 



T_ h _ J \v k Uk J \dk J ' \Vk I \ , l E k -vk 

k ' 7 \ \l 2E k 



The ground state \G) of this gapped system can be determined by the conditions j ks \G) = for s = ±1, which 
leads to the following form (not normalized) 



\G) = cxpS^—cldk + Y,— dlck\\G L )®\G R ) 
[k>o Uk k<o Vk J 



(16) 



with |Gl(_r)) the ground states of the gapless systems Hl and H R , respectively. The operators c k d k and d\c k creates 
particle-hole excitations in the massless theory. The physical meaning of this state can be understood better by 
redefining the quasi-particle creation operators in the massless theory as L\ k = c k and L 2k = c_ fc for k > 0, and 

similarly Ri k = —d_ k and i?2fe = d k for k < 0. In term of these operators, 



(17) 



with 

4 = + L\ 2 Rl k2 = - (c\dk + dl k c- k ) (18) 

creates a pair of quasi-particles in the two edges. Defining 

fe>0.s=l,2 



\ E ^(^f^)(4^-44+c_ feC t_ fe -^g 

fc>0,s=l,2 V 7 



(19) 



we have the identity 



Consequently \G) can be rewritten as 



\G) = e- H ' e w\^*l\e- H '\G L )®\G R ) 
I fc>o J 

= e- H ° exp { £ AU l G ^) ® l G «) = e ^ l G * 
U>o J 



(21) 



In the last line, we have used the fact that H e \Gl) ® \Gr) = 0. In the long-wavelength limit k — » 0, H e in Eq. (|X9 
is expanded as 

fc>0,s=l,2 s 

= ^ { - Hl+Hr) (22) 
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Consequently, Eq. (|21[) in the long wavelength limit correctly recovers the general relation ([6]) with the extrapolation 
length t = 1/2E S . 



